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We propose a class of curvaton models which we call passive curvaton. In this paper, two kinds 
of passive curvaton is considered. The first one is a pseudoscalar curvaton couples to a gauge field. 
Different from the inflaton case, the constraint from formation of primordial black holes (PBHs) 
is much weaker and large non-gaussianity (of the equiliteral type) can be produced. The second 
model is a dilaton-like scalar curvaton couples to a gauge field. We investigate the scale dependence 
of non-gaussianity in this model. In both models, the spectrum and non-Gaussianity are enhanced 
by the slow-roll parameter of the curvaton field. Other possible passive curvaton models are also 
mentioned. 

PACS numbers: 



I. INTRODUCTION 

Inflation has became very popular as a theory of the very early universe (for a textbook, see [l| and references 
therein). Unfortunately we still do not really know what causes inflation or what is the inflaton if inflation is driven 
by a scalar field. For single field inflation case, the universe is well discribed by the 'separate universe' approach 
[2] . Namely each Hubble patch can be regarded as a small Friedmann-Robcrson- Walker (FRW) universe and each of 
them 'looks the same' (adiabatic condition) in spite of a possible 'time shift' of its history 1 . The 'difference' of the 
number of e- folds of their evolution would corresponds to the curvature perturbation (so-called SN formalism 0, Q , 
see Q for a recent review). In this picture curvature perturbation is a constant (to any order of perturbation theory) 
on superhorizon scales due to the identity of each of the 'small universe'. Indeed, if each small universe evolves in 
the same way, intuitively an extra number of e-folds is not going to change with time. However, this picture can 
be changed if there exit another field, say a curvaton. Intuitively the independent fluctuation of the curvaton field 
will give each of the small universe an additional (and independent) contribution of energy density therefore they 
are not identical anymore (non-adiabatic) and the super-horizon curvature perturbation will evolve with time. If the 
condition is properly chosen, the curvature perturbation will increase from practically zero to the value we need to 
explain CMB temperature fluctuation and primordial density perturbation. This is what happens in the case of the 
curvaton scenario . Therefore initially the fluctuation of the curvaton corresponds to isocurvature perturbation 
(which means no curvature perturbation) at the beggining and finally the isocurvature perturbation is transformed 
to curvature perturbation. 

For single-field slow-roll inflation with inflaton field </>, the spectrum is given by 



where the slow-roll parameter is 



and Mp = 2.4 x 10 18 GeV is the reduced Planck scale. 
On the other hand, for a curvaton field a , 

pV 2 = r—=rr) a ^5a, (3) 
a a 

where a is the field value of curvaton at horizon exit and 0.01<r<lis roughly the ratio of the curvaton energy 
density to the total energy density of the universe when it decays. The lower bound of r is from the constraint of 
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1 This is an intuitive way of explanation but of course we actually has to be careful about the issue of gauge in general relativity. 
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non-Gaussianity of Cosmic Microwave Background (CMB) from WMAP The second equality in Eq. ((3]) holds for 
a quadratic potential of o with mass m and r\ a = m 2 /3H 2 is a slow-roll parameter of the curvaton. This form will 
become useful in the following sections. 



II. PASSIVE CURVATON 



Let us consider a curvaton field o couples to a certain quantum field \. The Lagrangian that is relevant to o is 
given by 



C=-g^d fl ad v a-V(a)+C I (a, X ), (4) 



where V is the curvaton potential, £/ is the interaction term, and the metric is 

ds 2 = dt 2 - a 2 (t)dx. 2 = a 2 (dr 2 - dx 2 ), (5) 
where r is the conformal time. The equations of motion for the mean fields arc then given by 

o + 3Ho + V = (6) 
do 

where p x is the energy density of \ particles, the dot and the prime denote differentiating with respect to t and a 
respectively. In Eq. ([6]), the right hand side of the equation is the back reaction of the interaction to the inflaton 
mean field. The back reaction arises due to copious production of x quanta during inflation. 
In addition, the fluctuations of o satisfy 

V 2 t , rrtlr Jdd 



So + 3HSo - \5o + V"So = S[^-). (7) 
a A \ do ) 

The homogeneous solution of this fluctuation equation gives rise to the standard primordial density fluctuations. Here 
we call them as active fluctuations and denote their power spectrum by P a . The right hand side of Eq. (|7|), which 
comes from the fluctuations of the backreaction, acts as a source for generating additional fluctuations of </>. The 
particular solution of Eq. (0 with the source term is referred as passive fluctuations and the power spectrum is denoted 
by Pp. Hence, from Eq. ([3]) the total power spectrum is given by the contributions from both active fluctuations and 
passive fluctuations as 

P = P a + P P = (-Y (\6o k \ 2 ), (8) 



where Sok is the Fourier mode of So and P a is given by 



r\2 / H x 



p *=u u ■ ^ 



For comparison, for an inflaton we would have (from Eq. ([I])) 

P=P a +P p = (S-) (\SM 2 ), (10) 



One benefit to consider the curvaton scenario is, different from the inflaton case, in Eq. ((7|), we do not need to worry 
about the effect of scalar metric perturbations because the energy density of the curvaton is subdominant at horizon 
exit and the fluctuation corresponds to isocurvature perturbation. 



III. PSEUDOSCALAR CURVATON 



For our first example, we consider the case in which the curvaton o is a pseudo Nambu-Goldstone boson (PNGB) 
with a typical potential after the shift symmetry is broken, 



V(o)=A i [l-cos(o/f)}, 



(11) 
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where A is a mass scale and / is the axion decay constant. We consider a coupling to a gauge field given by 

Ci = —vF^F^, (12) 

where = d^Ay — dyA^ is the field strength associated to some E/(l) gauge field A^ and F^ v is its dual. 
In order to calculate the spectrum, we substitute Eq. (p~2]) into Eq. ((?]) and obtain 



—^+2H ir -W 2 + a 2 m 2 



6<7(t, x) = a2 j [ E ■ B - ( E ■ B 



(13) 



where we have used conformal time t and H = da/ (adr). The homogeneous solution corresponds to active fluctuation 
and the particular solution corresponds to passive fluctuation. Working in Coulomb gauge, the circular polarization 
modes obeying [To[ 



dr 2 



k 2 ±^ 



A±(r,k) = 0, 



(14) 



where 



act 
2H] 1 



(15) 



and dot denotes differentiation with respect to t. There is a growth solution of fluctuations described by [T(| 



A+(r,k) 



1 / k 



1/4 



M-2y/^k/(aH) 



(16) 



This expression is valid for ^ > O(l). 



For the inflaton case , equation similar to Eq. ()T3|) has been solved [11|, ll2| 3 by using Eq. (|T6|) and the spectrum 



is (approximately) given by 



P = P a 1 + 7.5 x 10 



-5 



H 2 ( a\ 2 1 



Air 2 \2fJ e £ 



(17) 



where 



Pa 



ML 

2nd 



(18) 



is the active spectrum and 



We can relate Eqs. (JUJ) and (O via 



and write Eq. (fTTj) as 



4>a feaMp 



P, 



2Hf V 2 / 



H 2 f a \ ' 1 



a 4tt2 \ 2fJ i 



(19) 



(20) 



P = P a [ 1 + 7.5 x 10" 5 P a - 



(21) 



2 In this paper, the inflaton case refers to the case where \ couples to the inflaton instead of curvaton. 

3 The calculation is rather lenthy and there is no point to repeat it here. We refer the reader to the original references. 

4 We use the same notation for the inflaton case and curvaton case for comparison, but it should be clear through the context concerning 
which one we are refering to. 
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However, this result cannot be applied directly to curvaton. The reason is that the second equality of Eq. ((l~9|) no 
longer holds true because curvaton does not drive inflation. This fact will make a crucial difference between curvaton 
and inflaton case. 

The active spectrum of curvaton is given by Eqs. © and © as 

Vj£*- (22) 

The relation between curvaton active spectrum and £ is 

^ = 4^(27) e' (23) 

Note that this form is different from Eq. ([20)1 . The spectrum of the curvaton is given by substituting Eq. ([25)1 into 
Eq. (fTT]) with P a is given by Eq. ([22]). 

7.5 x 1CT 5 e 4 ^ 
p 

r 2 rfc ° £' 



P = Pa[i+ • -a-, Pa-ar). (24) 



The reason Eq. (ITTl) can be used here for the curvaton (with different P Q ) is because the field fluctuation satisfy the 
same equation. The difference is the way to 'convert' the spectrum of field fluctuation to the spectrum of curvature 
perturbation is different between inflaton and curvaton cases (see Eqs. |(5J) and (fTO)) ). 

For both inflation and curvaton case, the active fluctuatin should dominates the spectrum (becasue passive fluctua- 
tion is highly non-Gaussian) and CMB normalization should be imposed, namely, P a = 5 x I -5 Q for both inflaton 
and curvaton cases. Similarly we find the (equiliteral type) nonlinear parameter is given by 

/X'~4.4xl0 10 -^^. (25) 

As we can see from Eq. (|2~4")) , the passive fluctuation (corresponds to the second term) in the curvaton case is 
enhanced by the slow-roll parameter rj a and r. However, in general £ is expected to be smaller than the inflaton case 
because the Hubble parameter is determined by the inflaton potential instead of the curvaton potential therefore the 
enhancement is not necessarily large. But it does implies the effect of passive fluctuation can still be significant with 
a relatively smaller £ compared with the inflaton case. 

As pointed out by [l3[, for the inflaton case, primordial black holes may form. It can be seen from Eq. (flUl) that 
£ which appears in the exponent of passive fluctuation is proportional to y/e. However, for the active fluctuation of 
inflaton, as shown in Eq. (QJ, it is proportional to the inverse of y/e. In order to have slow-roll inflation at all, we 
need e< 1 during inflation. But inflation must come to an end therefore e is expected to increase toward the end 
of inflation. This implies although it is difficult to have PBH formation from active fluctuation, it would be possible 
for passive fluctuation. The constraints from PBH formation could be stronger than non-gaussianity which means 
non-gaussianity may be negligible. 

However, this is not the case for the curvaton becasue there is no direct relation between £ and e like in Eq. (|I9[) 
and unless a increases dramatically near the end of inflation which is not expected to happen, PBH will not form and 
large equiliteral shape non-Gaussianity can appear without constraint from PBH formation. 

IV. SCALAR CURVATON 

For a second example of passive curvaton, we consider 

Ci = -L^Lf^f^, (26) 
where I(<j) plays the role of a field dependent gauge coupling. Now the equation we need to solve is 

d 2 _a 



dr 2 dr 



2Utt- -V 2 + a 2 m 2 



<5<t(t, x) = a 2 ^j-[E 2 — B 2 }. (27) 



This is more difficult to solve than Eq. (fT5|) because the source term contains function of a. The inflaton case is 
considered in fl4j . However, because the inflaton dominates the energy density of the universe during horizon exit, 
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FIG. 1: Non-Gaussianity Jnl corresponds to number of e-folds N. The blue (upper) line corresponds to Ntot = 200 and purple 
(lower) line corresponds to Ntot = 100. 

the fluctuation of the inflaton causes fluctuation of the metric tensor therefore Eq. ([27]) have to be replaced by a more 
complicated version to take care of the metric perturbation. 

For the inflaton case, in order to simplify the calculation, we can choose the inflaton potential V(cf>) = M 2 (ft 2 /2 
and / oc exp(— (ft 2 ) to make I (x a 2 . For the curvaton case, we may assume both the inflaton potential and curvaton 
potential are quadratic, therefore we have 3H(ft = —M(f> and 3Ha = —ma and by deviding two equations to eliminate 
t we can show (ft oc a. Consequently we might also choose I oc exp(— a 2 ) to make I oc a 2 . This simplification also 
allows us to solve the vector modes. This exponential form (dilaton-like interaction) of / is typical in string theory 
and supergravity. For the inflaton case we have 



ht = L = 2 i L 

i lift (ft 



H 



[pl/2 



(28) 



where we have used I < 



Once again the last equality does not hold for the curvaton and we have 



I 

lb 



H 4nM P Pa 



1/2 



H 



rrja 



(29) 



From Eq. (f2"T)). we can see these difference makes the source term and the particular solution which corresponds to 
passive fluctuation different. 

The curvature perturbation in the inflaton case is given by fl4| 



P = P a [l + 192P a N 2 (N tot -N)] 



(30) 



For curvaton, we obtain 



P = P,. 



192 

l + - r - s P a N 2 (N tot -N) 



(31) 



and for the squeezed limit, non-gaussianity close to the local type is given by 
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Inl — 



0.7 



N 

60,' ( ^ 



N). 



(32) 



Again for the curvaton case, the passive fluctuation is enhanced by r\ a and r. However, for N ~ 60, the requirement 
that the first (standard) term in Eq. (f3~Tj) dominates implies that about Ntot < 600 for the inflaton case and N tot < 
600r 2 7y 2 r for the curvaton case. The constraint of total number of e-folds becomes more severe for the curvaton case. 
We plot the result in Fig. [TJ As we can see from the plot, the current constaint of non-Gaussianity /jvx < 100 
from WMAP @ requires N to t < 200 for ri] a ~ 1 and getting more severely for smaller rr\ a . The plot also suggests 
possibility to have a large running of Jnl- 

However, it can be easily checked that in this case PBHs is not generated near the end of inflation. Indeed when 
N ~ the second term in Eq. (|3"Tj) vanishes. 
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FIG. 2: Running non-Gaussianity at horizon exit corresponds to number of e-folds N. The blue (upper) line corresponds to 
Ntot = 200 and purple (lower) line corresponds to Ntot = 100. 

For the case of scalar curvaton, although we do not have PBHs formation after inflation, as the inflaton case in [l4| , 
large running of non-Gaussianity may be possible. We will investige this issue in more detail with recent observational 
constraints. 

If we allow the scale dependence of non-Gaussianity. f^L can be parameteriscd by using the running index n as 



f NL (k) cx W 



(33) 



or 



din I 



NL 



1 df. 



NL 



dink 



NL 



dN 



(34) 



By using Eqs. (|32|) and (f34|) . we can calculate n and the result is shown in Fig. [2] For N ~ 60 the result of 
running is quite predictive and does not depend on rr] a . We obtain roughly 0.01 < n < 0.1 and 30 < /jvx < 100 
for 100 < Ntot ^ 200. Although WMAP team assumes a scale independent Jnl, it is possible to analyse the data 
assuming varying Jnl- It is shown in that WMAP data allows running non-Gaussianity as large asn~ 0.1 5 with 
Jnl ~ 50 which is quite close to the result obtained here. The constraint is expected to be improved in near future 
for example by PLANCK satellite. 



OTHER POSSIBILITIES 



The idea of passive curvaton is not restricted to curvaton couples to gauge field. 

EMI, 



We can consider for example 



(35) 



where g is a coupling constant and <Ji is a constant field value. When cr rolls down to each trapping point at Cj , the Xi 
particles become instantaneously massless and are produced with a number density that increases with cr's velocity. 
As a dumps its kinetic energy into the Xi particles, it is slowed down and the produced Xi particles are diluted due 
to the inflationary expansion. However, due to the non-trivial time evolution of cr, we cannot apply Eq. ([3]) directly. 
It would be interesting to consider the detail of this model, but we will not persue it in the current work. 

Another possible direction is to consider models which involve only gravitational interaction without a direct non- 
gravitational Ci between curvaton and other quantum fields. For example we can consider the effect of quantum 
stress tensor fluctuations 120, Eli ■ 



To be more accurate, n = 0.3^^2 at 95% confidence. 
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VI. CONCLUSION AND DISCUSSION 

In this paper we consider two models of passive curvaton. The first one is a pseudoscalar curvaton and the second 
one is a dilaton-like scalar curvaton. We compare the models to the cases where the gauge field couples to an inflaton. 
We found in general the passive fluctuation is enhanced by the slow-roll parameter of the curvaton. In the first case, 
PBHs is not expected to form after inflation and non-Gaussianity can well be large. In the second case, we analyse 
the running of non-Gaussianity. 

Passive curvaton in principle can also be applied to the case of inflating curvaton [22| (see also [HI Hi!) The original 
inflating curvaton scenario predicts negligible non-Gaussianity because the curvaton dominates the energy density of 
the universe soon after curvaton start to inflate. However, the idea of passive curvaton can generate non-Gaissianity 
for inflating curvaton. 

For simplicity, we only considered curvaton field fluctuation to first order. However, in principle if we consider the 
fluctuation of curvaton field up to second order of field fluctuation, this together with passive flucutation can result 
in a non-Gaussianity of a complicated shape. However, the effects is subdominant and unlikely to be able to test in 
near future experiment. 
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